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FLOCKING AND NON-FLOCKING BEHAVIOR IN A 
STOCHASTIC CUCKER-SMALE SYSTEM 


TA VIET TONt, NGUYEN THI HOAI LINHt, ATSUSHI YAGI* *t 

Abstract. We first present a new stochastic version of the Cucker-Smale 
model of the emergent behavior in flocks in which the mutual communication 
between individuals is affected by random factor. Then, the existence and 
uniqueness of global solution to this system are verified. We show a result 
which agrees with natural fact that under the effect of large noise, there is no 
flocking. In contrast, if noise is small, then flocking may occur. Paper ends 
with some numerical examples. 


1. Introduction 


Flocking is a prevalent behavior of most population in natural world such as 
bacteria, birds, fishes. It is also widespread in some phenomena in physics, for ex¬ 
ample interacting oscillators. Recently, a number of articles proposed mathematical 
models for flocking behavior, to name a few [nisiiTiiiaiiiiiiiiiiTiiiH]. Vicsek and 
collaborators [TB] presented a model (for convenience, we call it Vicsek’s model) 
and then studied flocking behavior via computer simulations. Some theoretical re¬ 
sults on the convergence of that model can be found in [15] . Based on the Vicsek’s 
model, Cucker and Smale introduced a model for an A-particle system |4l[5]. Then, 
some mathematicians called it the Cucker-Smale system [HEmnidiKisKn]. we 
will recall here some features of this system. We consider motion of N particles 
in the space (d = 1,2,3,.. .). The position of the i-th particle is denoted by 
Xi = Xi{t) [i = 1,2,..., N). Its velocity is denoted by Vi = Vi{t) {i = 1,2,..., N). 
The Cucker-Smale system is as follows 



( 1 . 1 ) 


Here, the weights ~Xi \\) quantify the influence between i-th and j-th particles. 
This communication rate is a nonincreasing function '0: [0, oo) — >• [0,c») of the 
distances between particles. This function has various forms. In H m, iIj{s) = 
while in [III US], -ipis) = V’(s) = ^ or '0(s) = constant. For 

such functions, it is shown that when /3 < ^ the convergence of the velocities to 
a common velocity is guaranteed, while for /3 > ^ this convergence is guaranteed 
under some condition on the initial positions and velocities of particles. We call 
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them unconditional flocking and conditional flocking, respectively. In the latter 
case, the result on the non-flocking for two particles on a line is also stated [5]. 

We know that real systems are often exposed to influences that are incompletely 
understood. Therefore extending these deterministic models to ones that embrace 
more complicated variations is needed. A way of doing that is including stochastic 
influences or noise. Up to our knowledge, however, there are few such models which 
are studied theoretically [611I311I]. In , Cucker and Mordecki modified the model 
(HI) in by adding random noise to it 


\x', = Vi, l<i<N, 

Here communication rate "ip has the same form as that in (HD. Hi(t) is a three- 
dimensional Gausian centered, stationary stochastic process, that satisfies a d - 
dependence condition for some d > 0, i.e., two sets {i?i(s)|s < t} and {i?i(s)|s > 
t -I- (5} are independent for each t. And Hi{t) has trajectories and independent 
coordinates. The authors showed that a conditional nearly flocking occurs in 
finite time with a confidence which is similar to the conditional flocking of the 
deterministic system (HD- That is, for 1 / > 0 small enough, there exists a time To 
depending on i/ and initial values Ui(0) {1 < i < N) such that for every t G [0,To), 
Sij Ikj(^) ~ ^ ^ with a positive probability. This probability, however, is 

not one and the interval time for the occurring is not [0, oo). In other words, it is 
just nearly flocking, not flocking. In another approach, Ha and collaborators HU] 
studied the following two Cucker-Smale systems with the presence of white noise: 


dxi = Vidt, I < * < 

dvi = J2f=i fPiWxj - Xi\\){vj - v,)dt + y/Ddwi, 


and 

I dxi = Vidt, I < i < A, 

[dvi = A - Xi\\){vj - Vi)dt + D(vi - Ve)dwi, 

where is a constant vector. In |13] . the result on the flocking behavior of system 
(HD is exhibited when the communication rate is a constant. Furthermore, if the 
communication rate satisfies a lower bound condition, then the relative fluctuations 
of velocities around a mean velocity have a uniformly bounded variance in time. 
In [1], by giving another definition for flocking which is relative to almost surely 
convergence, the authors showed flocking of system (HD which covers the case of 
the communication weight employed by Cucker-Smale. 

We are interested in the communication rate term in system (HD- What hap¬ 
pens if this mutual communication is affected by random factor, i.e., ipilxj — 2 :^ 1 ) 
ip{\xj — Xi\) -I- white noise? This motivates us to present and study a Cucker-Smale 
system under the effect of a common white noise on the mutual communication. 
We will show that if the noise is large then there is no flocking. Even if the com¬ 
munication rate satisfies the unconditional flocking for (HD, this still holds true. 
This is different from the deterministic Cucker-Smale system (HD but is adaptive 
to real situations. For example, under a strongly random effect of strong winds or 
water currents, birds or fishes will separate and can not make a flock or school. In 
contrast, we will show that if the noise is small enough, then flocking occurs. For 
the convenience of calculation, we will use the Stratonovich stochastic differential 
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equations. Our model to study in this paper has the form 


dxi = Vidt, 


< l<i<N. 

[dvi = ° dwt, 

(1.5) 

Here cr > 0 is the strength of white noise and {wt,t > 0} is one-dimensional 
Brownian motion defined on a complete probability space with normal filtration 
(O, P). We assume that communication rate function ■!/):[0,oo) —>■ 

[0, oo) is locally Lipschitz continuous. 

The organization of the paper is as follows. In the next section, we prove the 
existence and uniqueness of global solution to (USD. In Section 3, non-flocking under 
the effect of large noise is shown. In contrast, in Section 4, we exhibit flocking under 
the effect of small noise. Finally, some numerical examples are presented in Section 

5. 


2. Existence and uniqueness of global solution 

In this section, we shall prove global existence of solution for the system CD). 

Theorem 2.1. For any given initial values (a;i(0), 11^(0)) £ (1 < * < system 

(irsi) has a unique and global solution. 

Proof. Since the functions on the right side of (lESl) are locally Lipschitz continuous 
on there is a unique solution {xi{t),Vi{t)) {1 < i < N) defined on an interval 
[0, r), where r < oo and if r < oo it is an explosion time [H H], i.e., 


Put 


T = sup{t: sup [||a:i(s)||-I-||ui(s)||] < oo}. 

se[0,t),l<i<A 


a:|| 


N 



N N 

. ^||u,||2, x=-Y,x,, 

\ 2=1 2=1 



2=1 


It follows from system (1131) that 


dx = vdt, 
dv = 0. 


Then almost surely x{t) = x{0)+v{0)t and v{t) = z;(0) for every t £ [0, r). Without 
loss of generality, we may assume that a;(0) = i;(0) = 0. Then, 

N N 

'^Xi{t) ='^Vi{t) = 0 a.s. (2.1) 

2=1 2=1 


^ \\v,{t)-vm" = 2N\\v{t)\\^ 

ij'=l 


( 2 . 2 ) 
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It follows from the second equation of (II.hp and from the chain rule of Stratonovich 
stochastic differential equation that 


N 


d\\v\f ^11^*11^ 

N 

=2^ < Vi,dvi > 


2=1 

N N 

—2 ^ V’dki — ^dl) ^ — Vi > dt 2a ^ < Vi,Vj — Vi > odwt- 

ij'=l ij'=l 

(2.3) 


We have 


N 


N 


N 


< Vi,Vj - Vi > = <Vt- Vj,Vj - Vi > + '^ < Vj,Vj - Vi > 


i,j=i 


which induces 


* j=i 
N 


hi=i 


N 


= - X! 11^* “ “ X! <Vk,V^- Vk >, 




2,fc=l 


N 

2E 


AT 


= -E 


Ui Uj I 


ij = l 


= -2Af||u|p (seedU)). 


(2.4) 


Furthermore, it follows from 

N 


N 






i,t=l 


i,j=i 

N 


+ E > 

*0 = 1 


that 


N N 

2 E -a;*||) < >= - E -2;i||)||ui -Ujlp. (2.5) 

i,i=i hi=i 

Thus, by (I2.3D - (I2.5I) . we obtain 

N 

d\\v\\‘^ = - tp{\\xj - a;j||)||uj - VjW'^dt - 2Na\\v\\'^ o dwt- 
i,j=i 

Or, equivalently, in the ltd form: 

N 


d\\vr = 


- E V'dki -a;^ll)ll^'^ -olP + 2A^V2||u| 

*,i=i 


dt — 2iV(T||u|pdw(. 


( 2 . 6 ) 
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Hence, by using the comparison theorem [14], it follows from (I2.6p that for every 
t>0, lk(^)lP ^ ^(.t) a.s., where V{t) satisfies the following equation 

idV = 2N^a^Vdt - 2NaVdwt, 

\F(0) = |K0)|p. 

This linear equation has a unique global solution V{t) = . Thus for 

every t G [0, r) 

IK0ll<IK0)||e-^'^- a.s. (2.7) 

Then, from the first equation of (1I3|) we have almost surely 

N N 

< Xi,Vi > dt 

i=l i=l 

< 2||a:|| \\v\\dt < 2||w(0)||e“'^‘^“*||x||dt. 


By the comparison theorem, we obtain ||a;(t)|p < X{t) for all t > 0, where X{t) 
satisfies the following equation 

jdX = 2||u(0)||e-"‘"“*\/Xdt, 

1x(o) = |Ko)|p>o. 


Since X{t) = [||a;(0)|| + ||u(0)||/q e then for every t G [0,r) 


x(t)||<||x(0)|| + |K0)|| 



e-N<y^.ds 


a.s. (2.8) 


From (EJj), (EH) and the definition of r, we see that t = oo a.s. It means that the 
solution to (EH is unique and global. □ 


3. Non-flocking under large noise 

In this section, we will show that under the effect of large noise on particles, 
i.e., a is large, then the system EH does not flock. First, we give a definition for 
flocking. Then non-flocking theorem is shown. 

Definition 3.1. The state of particles Vi{t)) {1 < i < N) in system (11.51) has 

a time-asymptotic flocking if, for 1 < ifl < N, the velocity alignment and group 
forming in the following senses, respectively, are satisfied 

1. limt_>oo = 0. 

2- supo<t<^ E| 11 < oo. 

Put a = supg>Q i/'(s). Throughout this section, we assume that the communica¬ 
tion rate satisfies an upper bound condition, i.e., a < oo. Note that the communi¬ 
cation rate in Cucker-Smale system satisfies this condition. 

Theorem 3.2 (Non-flocking theorem). If a > then the particles do not flock. 

Proof. As we see in the proof of Theorem 12.11 that without loss of generality, we 
can assume that iE(0) = h(0) = 0. Then (12.11) and (12.21) hold true. It follows from 
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(I2.2tl and dill) that 

N 

d\\v\\'^ = ^ - tp{\\xj - a:i||)]||i;i - VjW'^dt - 2Na\\v\\'^dwt 

ij=l 
N 

> ^ (-/Vct^ — a)||z;j — VjW'^dt — 2Na\\v\\'^dwt 
ij=l 

=2N{Na‘^ — — 2Na\\v\\‘^dwt. 

By using the comparison theorem, for every t > 0 we have ||i’(t)|P > 
where Vi(t) is a unique solution of the following equation 

f dVi = 2N{Na^ - a)Vidt - 2NaVidwt, 

\l/i(0) = ||u(0)||2 >0. 

Since dEVi = 2N{Na‘^ — a)EVidt, we have EVi(t) = —>■ oo as 

t —>■ oo. Therefore, E||u(t)|p —>■ oo as t —>■ oo. It means that the particles do not 
flock. □ 


4. Flocking under small noise 

In this section, we consider the system m under the effect of small noise, i.e., 
a is small. We will show that flocking takes place under a lower bound condition 
of the communication rate. 


Theorem 4.1 (Flocking theorem). Assume that there exists ip* > 0 such that 

inf ip{s) > Ip*. 

s>0 


If a < 
i < N). 


^ then particles flock under any initial values {xi{0), Vi{0)) 


p2d 


(1 < 


Proof. Similar to Theorem 13.21 we can assume that a:(0) = u(0) = 0. Then (12.11) 
and (|2.2I) hold true. By (12.21) and (12.61) . we have 


N 


^ 'f’*\\xi — VjW'^ + 2iV^(T^||u|p 
id 

= (^Na'^ — i/)*) ||u|pdt — 2A^cr||u|p(iu't. 


dt — 2N a\\v\'\^ dwt 


By using the comparison theorem, it follows that ||u(t)|p < V 2 {t) a.s., where V 2 {f) 
satisfies the following equation 


dV2 = [Na^ — ip*)V2dt — 2NaV2dwt, 

V 2 { 0 ) = \\vm\f 


Since fiEF 2 = (fVcr^ — ip*'jlE,V 2 dt and a < \J~^^ we have 

]E|k(i)|P < EV 2 {t) = ^ 0 as t ^ oo. 

On the other hand, the linear equation (HU) has an explicit solution 

V2{t) = \\y{^Q)fe^(N-2N^)o^-r]t-2Naw,^ 
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Thus, 

\\v{t)\\< (4.2) 

In the prove of Theorem 12.11 we showed that d||a:|p < 2||a;||||u||(it. Then by using 
the comparison theorem and (14.21) , it is easy to obtain that 

||a;(<)|| < lk(0)|| + ||u(0)|| [ 

Jo 

Thus, 


E||a:(t)|| <||x(0)|| + ||i;(0)|| 
= ||x(0)|| + |K0)|| 
= ||x(0)|| + |K0)|| 
= ||x(0)|| + |K0)|| 


^^i[(Ar_2A=),T=-V>*]SgiJV^rT=.^^ 




{N — N‘^)a^ — Ip* . 


j[(Ar_Ar2)^2_^.]j 


1 


Therefore, supj>gE||x(<)|| < oo. From the above results, we conclude that the 
particles flock. □ 


Remark 4.2. The lower bound condition has been used for the system (ESI) by 
Ha et al. m- The authors showed that when this condition holds true, the term 
Si jE’llui — WjlP is uniformly bounded in t. Furthermore, when '0 = 1, flocking of 
particles takes place. However, when 0 satisfies the lower bound condition but is 
not a constant, the flocking result has not yet obtained in m- 

Remark 4.3. The above flocking result is relative to the expectation (average) of 
solution. We can obtain a result on flocking of almost surely convergence that is 
valid for individual trajectory. To see that, firstly, we change the first condition in 
Definition 13.II by limt^oo H'l'i ~ '^ill =0 a.s. Now, it follows from (14.21) that 

\\vm < ||z;(0)||e*[5(^-2^^V^-b*-^-^]. (4.3) 

By the strong law of large numbers for Brownian motion [9], ?>0ast —>-00 a.s. 

Thus limt^.oo ||x(t)|| = 0. It then follows from (12.21) that for 1 < f < iV, 

lim llvi —i;i||=0 a.s. (4.4) 

t-yoo " J" 

It means that flocking occurs. Furthermore, it is easy to see from gSD that (I4a 
still holds true for the case cr > 0, '0* = 0. This means that flocking occurs without 
the assumptions on the positivity of '0* and upper boundedness of a. Consequently, 
unconditional flocking occurs not only for the communication rate 0 in |13) but for 
those in IlISJIS]. 


5. Numerical examples 

In this section, we present some results of numerical simulations based on Euler’s 
method for the system (E^. First, we give examples for non-flocking, second, 
examples for flocking. 
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5.1. Non-flocking. Let us first observe an example that shows that the system 
(USD does not flock. 

Set a = 1, cr = 0.3 and i/>(s) = compute 100 trajectories of solution 

of (11.511 in with TV = 50 and initial values (a;i(0),'(;i(0)) are randomly generated 
in [0,0.1]^. Figure 1 illustrates behavior of function f{t) = ~ 'Cj(^)ll^ 

up to time T = 0.2. We see that values of function f{t) are very large even at small 
time t. 



Figure 1. Non-flocking in 2-dim. space 


5.2. Flocking. Let us next observe another example showing flocking. Set a = 
l,tT = 0.05 and '0{s) = 1. Initial values (a:i(0),Wi(0)) are generated randomly in 
[0,1]®. Figure 2 shows a flock of 50 particles at T = 0,0.02,0.5,1. Each vector 
shows position and direction of motion of each particle. The lengths of vectors 
represent the magnitudes of velocity vectors of particles. 

6. Conclusion 

This paper presented a new stochastic version of the Cucker-Smale model in 
which the mutual communication between individuals is affected by the common 
random factor. This model can find its applications in real world because almost 
every real phenomena are subject to environmental noises. Under the upper bound 
condition of communication rate, we show that if the noise is large then the particles 
can not flock. This is consistent with some situations in nature, for example, under 
a strongly random effect, there is no flock. In contrast, under the lower bound 
condition of communication rate and small noise, flocking occurs. The paper ends 
with some numerical examples of both flocking and non-flocking. 

Using the sense of flocking stated in Remark 14.31 unconditional flocking takes 
place. However, a remaining interesting problem for the stochastic Cucker-Smale 
model is to obtain flocking results in the sense of Definition 13.11 By numerical 
simulations, we predict that flocking in this case does not only depend on small 
noise but on some certain domain of initial values (a:i(0),?;i(0)). This problem is 
left for our future research. 
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T=0 



T=0.5 



T=0.02 




Figure 2. Flocking in 3-dim. space 
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